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The theory of bigravity offers one of the simplest possibilities to describe a massive graviton 
while having self-accelerating cosmological solutions without a cosmological constant. However, it 
has been shown recently that bigravity is affected by early-time fast growing modes on the tensor 
sector. Here we argue that we can only trust the linear analysis up to when perturbations are in 
the linear regime and use a cut-off to stop the growing of the metric perturbations. This analysis, 
although more consistent, still leads to growing tensor modes that are unacceptably large for the 
theory to be compatible with measurements of the cosmic microwave background (CMB), both in 
temperature and polarization spectra. In order to suppress the growing modes and make the model 
compatible with CMB spectra, we find it necessary to either fine-tune the initial conditions, modify 
the theory or set the cut-off for the tensor perturbations of the second metric much lower than unity. 

Initial conditions such that the growing mode is sufficiently suppresed can be achieved in scenarios 
in which inflation ends at the GeV scale. 


I. INTRODUCTION 

Evidence from an increasing number of cosmological observables favours an accelerating universe at late times m- 
This era of accelerated expansion may be due to novel gravitational physics, which will be tested by ongoing and 
future experiments m■ This possibility has triggered vigorous interest in alternative theories of gravity ®mm- 
Any modification of gravity requires new degrees of freedom (dof). Since the theory of a massless graviton is unique, 
new dofs are often gained by adding new fields. The simplest possibility is the addition of a scalar field, typically 
resulting in theories belonging to the Horndeski class mm] or beyond ii. 

Formulating a theory of massive gravity has been a long standing problem in theoretical physics due to the difficulties 
to incorporate the right degrees of freedom. The linear Fierz-Pauli theory had been developed long time ago [T7], but 
until recently all non-linear completions introduced the so called Bouleware-Deser (BD) ghost [TS], an extra dof that 
makes the theory not viable. Despite the difficulties, a class of healthy theories has been recently identified m in 
which a specific choice of the potential terms makes the theory ghost-free [20]. All these theories of massive gravity 
describe an interaction of two tensor fields in which the second one, the so called reference metric , is fixed. While 
massive gravity only allows static solutions on homogeneous backgrounds EH a bimetric theory with a dynamical 
reference metric does not introduce the BD ghost and describes dynamical cosmologies |22W21] (see also the reviews 
[25l 1 26] ). Cosmological solutions in these bimetric theories often allow for self-acceleration without the introduction 
of a cosmological constant m and were successfully compared to observations at background level |27fI2T7| . 

Many bigravity theories are however affected by gradient instabilities in their scalar sector, as has been shown 
by studies of the linear perturbations |30f!52] (see Refs. mm for derivations of the equations Refs. |.‘Vll - fiT%] for 
discussion of their dynamics). Stable evolution can be achieved only in a two parameters class of models known as 
Infinite-Branch Bigravity (IBB) [3F. In IBB, the reference metric (in keeping with common usage, we keep referring to 
the second metric as reference metric even if in reality is dynamical; we also use the notation /-metric) is contracting 
during the radiation and most of the matter era, until it undergoes a bounce at low redshift and begins to expand, 
coinciding with the onset of accelerated expansion in the physical metric without the need for a cosmological constant. 
The early time contraction of the reference metric makes tensor perturbations grow with time in IBB theories, as it 
was first shown in Ref. 136]; ;3jL (see also [3T[ 00] for modified tensor perturbation equations). This growing mode 
couples to the physical metric and severely modifies its dynamics, leading to observable consequences. 

In this paper we will investigate the effects of these large tensor perturbations on the Cosmic Microwave Background 
(CMB) and possible mechanisms to make the theory compatible with current observations. The perturbations in the 
reference metric grow very fast and rapidly become non-linear. At this point we will assume that tensor perturbations 
stabilize, modeling this effect by introducing a cut-off in the perturbations of the reference metric. Despite this 
treatment, the tensor growing mode significantly affects the evolution of the physical metric, and the consequences 
can be seen as an enhancement of both temperature and polarization spectra on low multipoles. These effects cannot 
be sufficiently reduced by varying the bigravity or other cosmological parameters: making the theory viable requires 
either fine tuning of the initial conditions, lowering the cut-off or modifying the theory. As it will be shown below, 
sufficient suppression of the growing mode can be achieved by an inflationary mechanism that produces Hubble-scale 
tensor perturbations at an energy scale of order few GeV. 
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II. BIGRAVITY 


We start with the action of the form [22] 

s = —f J d 4 x^gR(g) - J d 4 x^fR(f) ( 1 ) 

4 

+ m 2 M 2 f d 4 x^g ^ /3 n e n {X) + f d 4 x^gC m 

^ n =0 ^ 

where e„(A) are the elementary symmetric polynomials of the eigenvalues of the matrices A“ = M g and 

Mf are the Planck masses for g^ v and f gl/ , respectively, m is the mass scale of the graviton, /3 n are arbitrary constants 
and C m = C m {g ,V0 is the matter Lagrangian. Throughout the paper we will use a mostly plus metric signature 
convention and natural units in which the speed of light c is set to one. 

Here g^ is the standard metric coupled to matter fields in the C m Lagrangian, while / M „ is an additional dynamical 
tensor field. In the following we express masses in units of the Planck mass M g and the mass parameter m 2 will 
be absorbed into the parameters f3 n . Varying the action with respect to g one obtains the following equations of 
motion: 


G 


fllS 


-^(-1 )"m 2 /3 n 


71=0 


9»\Y$i) v (X) + gvxYfaftiX) 


= T, 
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( 2 ) 


where G is Einstein’s tensor, and the expressions PA (A) are defined as 


*(o) = G (3) 

Y (i) = X-I[X], (4) 

y (2) = A' 2 - X[X\ + \l ([A] 2 - [A' 2 ]) (5) 

y (3) = A 3 - A 2 [A] + \x ([A] 2 - [A 2 ])) 

-i/([A] 3 -3[A][A 2 ]+2[A 3 ]) (6) 


where / is the identity matrix and [...] is the trace operator. 
Varying the action with respect to f gl , we get 
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f,xY^ )u ( A" 1 ) + f vX Y^{ A” 1 ) 


= 0 


(7) 


where the overbar indicates f ^ curvatures. Notice that /do acts as a pure cosmological constant, which is however not 
required to satisfy current observations. Finally, the rescaling / —► Mj 2 f , /3 n — > MJ/3 n allows us to assume Mf = 1 
in the following (see |H]). Additionally, from now on we absorb the graviton mass rn into the constants /3*. 

We assume now a cosmological spatially flat FRW metric: 

ds 2 = a 2 (r) (—dr 2 + dxidx' 1 ') (8) 

where r represents the conformal time and a dot will represent the derivative with respect to it. The second metric 
is chosen also in a spatially FRW form 


ds) = - 


b(r) 2 /TL 2 (t) dr 2 + b(r) 2 dxidx 1 


(9) 


where H = a/a is the conformal Hubble function and b(r) is the ‘scale’ factor associated with the second metric /. 
This form of the metric ensures that the equations satisfy the Bianchi constraints (see e.g. [21]). 

The background equations for the two metrics have been obtained and discussed at length in several papers EH 
ESI® 02]. Here we summarize the main properties in the notation of (35] . Defining r(r) = 6(r)/a(r) as the ratio of 
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the two scale factors, the background equations can be conveniently written as a first order system of two equations 
for r(t) and 77: 


2 WW. + T-L 2 = a 2 (Bo + B 2 r' - w to tPtot), 
/ 3vBi^li 0 i(\ T w to t ) 


Pi — 3/?3r 2 — 2/?4r 3 + 3 B 2 r 2 ’ 


( 10 ) 

( 11 ) 


where the prime denotes derivative with respect to N = log a [22, 2D], w to t denotes the equation of state (EOS) 
corresponding to the total density parameter fl to t and the functions Bq(t), -Bi(t), B 2 (t ) are related to the Pi and r(r) 
as follows: 


B 0 (t) = p 0 + 3/3i r + 3/3 2 r 2 + far 3 , 
B\(t) = Pi + 3/3 2 r + 3 P 3 r 2 + ^r 3 , 
B 2 (t) = pi + 2p 2 r + p 3 r 2 . 


( 12 ) 

(13) 

(14) 


For simplicity the time dependence of So, 1,2 will be understood from now on. The Friedmann equation (i.e. the (0, 0) 
component of Eq. §) gives 

3H 2 =a 2 (p tot + B 0 ), (15) 

and by combining with the (0, 0) component of Eq. 0 we obtain a useful relation between 77 and the ratio r(r): 

a 2 B l 


= 


3 r 


(16) 


Finally, the combination of the last two expressions for % provide flt. 0 t( T ) = 1 — |q r ( r ) which can be inserted in Eq. 
(Ill to produce a closed differential equation for r(r) alone. 

The behavior of the background solutions depends on the choice of the pi constants and on the initial value of r. 
We denote solutions with the same Pi but different initial conditions as branches of the same theory. In Ref. [271 
it was shown that for each choice of Pi only two branches exist that agree with a standard cosmological early time 
evolution (like a matter dominated era at early times) and allow for physical solutions (e.g. p, 77 > 0). In the first 
branch, r evolves from r = 0 to a de Sitter point at a finite value r c > 0. These branches, however, suffer from scalar 
instabilities EH]. Only choosing P 2 = P 3 = 0 and the second type of branches in which r evolves from r —> 00 in the 
asymptotic past towards a de Sitter point at a constant r c > 0, allows for stable scalar perturbations and is compatible 
with background data [ 25 ): we dubbed this case infinite-branch bigravity (IBB). Note that even though an additional 
non-vanishing effective cosmological constant Po is viable, we assume Po = 0 since it would not affect the early-time 
evolution and is not required in order to fit observational data (see [25] 22]). From now on we restrict ourselves to 
IBB, in which only pi and P 4 are non-zero. This choice avoids introducing an explicit cosmological constant, which 
would make the entire bigravity model somewhat less appealing. 

As shown in [35], IBB models have to satisfy 0 < P 4 < 2p\ in order to get an initial value of r on the infinite branch. 
In particular, it was found that the best fit model occurs for pi = 0.48 and P 4 = 0.94: from now on we refer to this 
choice as the reference IBB model. We then have 


n 2 = 


a 2 (Pi + p4r 3 ) 

3 r 


(17) 


Here we derive the early time behaviour of the background evolution for later use. (corresponding to early time in 
IBB), Eq. (11) for IBB reduces to 


r ' ~ + Wto ^ r 

so that for wtot = const (i.e., in radiation or matter dominated epochs) one has 

r ~ 0 "3(i+™*„«)/2 5 

and the f^ scale factor b(r ) = r(r)a(r) goes as: 


(18) 


(19) 
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-(l+3w t ot)/2 


( 20 ) 
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The scale factor b(r) therefore contracts instead of expanding as long as Wtot > —1/3. Moreover, in the same 
approximation, 


n 2 


a 2 ftr 2 

3 


( 21 ) 


It is useful to derive an approximated estimate for the 6(r) bounce epoch. The bounce occurs b' = 0 <t=> r' = —r b - 

Tb 

If we assume a bounce after the radiation epoch, then the ratio at the bounce has to satisfy 

4ft - 6ft r 2 b + ftr? = 0. (22) 

From comparisons with observational data, we know that the best fit is close to ft ~ 2ft which leads to 

r b ~ 1 + yft . (23) 


It is also useful, for future purposes, to take note of the approximate observational relation between ft and 

n2 ^ (1 ^t.otO) 


ft = , 


D)-ft 2 


3ft 


(1 — ^toto) 


(^toid ~ 1)' 


+ <ftft 


(24) 


This relation approximately corresponds to the degenerate line between ft and ft for a flat universe, when fitting 
data sets such as supernovae, CMB and baryonic acoustic oscillation data [33]. Moreover, we require ft < 0.5, to 


ensure that the solution of Eq. (161 at present time lies on the infinite branch. 


III. TENSOR PERTURBATIONS 


As we are interested on the effect of bigravity on gravitational waves, we now proceed with writing the tensor 
perturbation equations ED ESI ESI- For the perturbed part of the metrics we adopt the transverse-traceless (TT) 
gauge, i.e. we select a transverse wave propagating along the z direction. Then, the tensor metric perturbations are 
given by: 


/ ft(+) ft(x) 0 \ 

hg(ij) = I ft(x) — ft?(+) 0 

\ 0 0 0 / 


(25) 


and similarly for the tensor modes of the / metric. We then obtain the following equations for both components 
(suppressing the subscripts +, x) 

h'n + 7 nh' n + (m 2 + (? n u~ 2 k 2 ) h n = q n h m , (26) 

where the indices n ^ m refer to ^-metric and /-metric, respectively; we have then 


. w 

^ = 2 + ^-’ ft 

m 2 = l-L~ 2 a 2 Br, m 2 

2 i 2 

ft = ft 

q g = T-L~ 2 a 2 Br , q f 


2r 2 + 3r' 2 + r (4 r' - r") 
r (r' + r) 

(r' + r) 2 n . 

H 2 r 2 

(■ r' + r) 2 


( r ' + r ) .ftp 


n ’ 


(27) 

(28) 

(29) 

(30) 


and where: 


B = pi + ftr 2 + r (2ft + ftr 7 ) + ftr 7 . 


(31) 


These equations are equivalent to the ones in Refs. ED ESI . In IBB (i.e. for ft = ft = ft = 0), B is simply given by 
ft. The coefficients (27>-(301 for the two tensor equations are plotted in Fig{l]as a function of redshift, for the choice 
ft = 0.48 and ft = 0.94. For this reference model, when considering both matter and radiation, the bounce happens 
at a redshift z b — 0.9, with a corresponding r b — 2.8. 
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FIG. 1: Coefficients of the tensor equations for IBB ( 27p0 l. Solid/dashed lines indicate positive/negative values and the vertical 
dotted line marks the bounce of the reference metric. Note that the friction term in the /-metric is negative at early times, 
when the /-metric is contracting. Note also that the coupling and effective mass terms are equal, and the ones corresponding 
to the / metric are very suppressed at early times. The bottom-right panel shows the evolution of the two scale factors. 


Let us anticipate here an important feature of these equations. As it will be shown below (see also Ref. |39]1. the 
equation for hf is unstable at early times since its friction term is negative as long as the scale factor b(t) is collapsing 
instead of expanding. The fast growth of hf will then drive a fast growth of h g as well, through the coupling term. 
However, in the limit r —> oo, the coupling coefficient q g becomes 


a 2 rpi 
qg = H 2 


r— >-oo v Q Pi 

y o - 


Pa r 


(32) 


and is therefore relatively small for large r. For the reference IBB model we have Pi/P a ~ 0.5; more in general, 
according to Eq. (241, Pi/P a ~ l/3/?i to within factors of order unity, and therefore since pi < 0.5, we have the lower 
bound 


% 


1 1 

- > - ■ 

Pir r 


(33) 


At recombination, for instance, we have r ss 10 4 in IBB so that one needs a hf roughly 10 4 times bigger than h 


before the coupling term q g hf becomes comparable to the h g terms in Eq. (26) and it starts driving the evolution 
of h g . This means that, in principle, a growing mode in hf will take some time before affecting h g . Whether this is 
enough to spoil the physical metric, is what we are going to test below. 
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In the following subsections we discuss more in detail the time behavior of h g ,hf during the inflationary, radiation 
and matter eras. 


A. Inflation 


During a de Sitter epoch in which H = const, one has %' ~ H and from Eq. (Ill: 

r ~ const . 


The tensor equations (261 then reduce to: 


h'q + 3 b! +hg (--2 + 


k 2 a 2 far 


-g ' ^9\ U 2 ' U 2 


) = 


'*5 + 3 7 + + ^) 


a 2 far 

a 2 13,, 

rU 2 9 


(34) 


(35) 

(36) 


We can now assume a 2 far <C H 2 during inflation (ie pi n f p m g) so h g behaves as in GR. The same is true for hf 
since a 2 fa/rhL 2 ~ {fa/fa )r~ 3 -C 1 . Since the inflationary equations are the standard ones, we expect the initial 
conditions to be unchanged and to apply equally well to h g and hf. 


B. Radiation and Matter Dominated Era 


In the early time, we can approximate the ratio of scale factors as r' = — | (1 + wtat) r which is solved by 

r = Aa~i {l+WM) , 


(37) 


where A is a suitable normalization constant of order unity. Furthermore we approximate 7~L 2 — I faa 2 r 2 . If the initial 
conditions for h g and hf are similar, then the source terms (301 are negligible at early times, i.e. small a, and the 
equations decouple. Furthermore, we find 


7 9 

~ f (1 - Wtot) , 

IS - ~2 ( 3wt °* + X ) > 

(38) 

m 2 

~ § (u> tot + l) 

- A/3 4 “ 

m 2 _ 3/3i(3 w tot + !) %(w tot +i) 

f ~ 2A 3 fa 

(39) 

4 

~i, 

2 (3 Wtot + l) 2 

c f- 4 

(40) 


Neglecting the mass term m 2 at early times, the tensor evolution for is described by 

h';-faw tol+ i)h ' f+ h < k ^ +1)2 =o. 


(41) 


At large scales the last term is negligible and one finds a growth of hf as a 3 ( 3u, t°t+ 1 )/ 2 . Thus, when radiation dominates, 
hf increases very fast as a 3 . Clearly, if one starts with h/ = 0 then this growing mode is initially absent and it takes 
some time before it becomes visible. The evolution of h g has instead a constant mode h ~ const until hf is large 
enough to source the growth of h g , cf Eq. (321. 


The early time approximation that leads to Eq (411 turns out to be a very good approximation also in the matter 
domination. In this regime hf increases as a 3 / 2 for super-horizon modes. When the coupling term becomes important, 
h g is driven by hf and acquires the same trend. Finally, when MDE ends and the system approaches a de Sitter 
behavior, the perturbations begin to decay. 

For sub-horizon scales the behavior is influenced by the hf time-dependent sound speed. An asymptotic form for 
large k can however be found. In this regime we can neglect the mass and the coupling terms, and the h g ,hf equations 
during either RDE or MDE have the general form 


h" + + fak 2 a v h n = 0 , 


(42) 
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FIG. 2: Evolution of tensor perturbations at different scales (top panels) and scale dependence of tensor perturbations at 
z = 0 (bottom left panel). All scales include a cut-off when non-linearity is reached; h g does not reach the cut-off for any 
scale and redshift. The initial conditions have been chosen so that hf = 1 is initially non-linear (solid lines) or starts at a 
small value (dashed lines) and becomes non-linear at latter times. The latter choice corresponds to h'f^ in f « h/(in) = 10 _2 ° at 
a = 1CP 10 (see section V Al. The two modes shown in the top panel correspond to the vertical dash-dotted lines in the bottom 
left panel. For reference, we recall the standard CM B ph oton visibility function (bottom right panel), whose peaks correspond 
to recombination and reionization epochs (see section IVI. The bounce of the reference metric has been indicated with a dotted 
vertical line. 


where the index n stands for g, f and r) = 1 + 3 Wtot and /3 n is an irrelevant constant. The general solution can be 
easily written in terms of the Bessel functions but here we need only the asymptotic behavior for large k or late times, 
which is 


1 \ 2 ' 4 / 


(43) 


times fast oscillations. We see then that for sub-horizon modes hf grows as a 1 in RDE, as a 1 / 2 in MDE and a final 
decay as a -1 when approaching the future deSitter phase, while h g decays as in the standard case as a -1 in all eras 
(before being driven to growth by the coupling to hf). For very large wavenumbers the coupling and the mass terms 
are ineffective at all times and the h g equation reduces to the standard case. This implies that there is no large effect 
to be expected for the directly detectable range of gravitational waves, which is around 0.1Hz or k « 10 14 Mpc -1 (see 
e.g. [33j), although a precise calculation is beyond the scope of this paper. 



































C. Inflationary initial conditions 


We can now use the a 3 growth mode during radiation to estimate the order of magnitude effect of the tensor modes 
at recombination (a more precise estimation will be obtained numerically in sections IV Vj. Inflation ends at some 
energy scale that can vary from 10 15 GeV to few MeV depending on the model. The upper limit comes from the 
bounds on the amplitude of tensor perturbations, indicating that the energy scale of inflation is at least that of Grand 
Unified Theories when observable modes are produced. The lower bound is inferred from the need of a radiation 
dominated universe in thermal equilibrium during big bang nucleosynthesis. These values are reached when the scale 
factor was a,i n f ~ 10 -9 at the latest. Since super-horizon tensor modes grow as a 3 during radiation domination, in 
the most favourable case of inflation ending just before big bang nucleosynthesis one would obtain an enhancement 
until recombination a rec ~ 10 -3 of hf ( lec ) ~ 10 18 h/( e ) roughly, where the subscript e denotes the end of inflation. If 
hf( e ) has the standard value approximately equal to H e /Tp « T^/Tp during inflation, where Tp ss 2.4 • 10 18 GeV 
is the reduced Planck temperature/energy and T e is the inflationary energy scale (here for simplicity assumed to be 
similar to the energy at the end of inflation), then the value at recombination of hf for a wave that reenters horizon 
at recombination or larger is roughly 


hf(rec) 




(44) 


Shorter waves reenter before and therefore grow less. If this value has to be compatible with the level of fluctuations 
in the CMB polarization spectra, then it should be lower than about one tenth of the temperature fluctuations; we 
take conservatively the level 1(V 7 . The same value should be taken for h'f since inflation excites both the tensor mode 
and its momentum conjugate. However we do not detect directly hf but rather the ^-metric mode h g which is coupled 
to matter, so as already noticed one can have a value of hf larger than h g by a factor of q ~ 4 w 10 4 at recombination. 
Putting therefore hf( rec ) < 10“ 3 , we obtain an upper limit to the temperature at the end of inflation T e ss 10 GeV. 
Since tensor modes impact CMB also at reionization, when the coupling term q g is closer to unity, this limit should 
be lowered to roughly 


T e « 1 GeV (45) 

(a similar limit has been obtained also in (23). It might be interesting to remark that the super horizon growing mode 
breaks the standard link between tensor modes and inflationary scale due to the presence of the coupling: now in 
principle one can have observable tensor modes even in low scale inflation. 

Any inflationary model with higher energy scale will generate excessive power on the tensor modes unless the 
inflationary initial conditions are suppressed with respect to the standard value or their growth is reduced. Taken at 
face value, this shows that the a 3 growing mode can be reconciled with observations only in the rather extreme scenarios 
of very low-energy inflation, as e.g. in the models discussed in Ref. (mas). Fixing the initial conditions to a more 
conservative era for the end of inflation, e.g. T ~ 10 3 GeV, would produce the huge value /i/( re c) ~ 10 39 /i/( e ) ~ 10 7 . 

Barring the case of very low-energy inflation, then, the IBB model is at odds with CMB observations. In the rest 
of the paper we will explore more or less contrived ways to overcome this difficulty. 


D. The non-linear cut-off 

The tensor perturbations in the reference metric grow so fast that they will eventually become non-linear. At 
this point, the perturbative treatment followed so far breaks down and one has to take into account higher order 
corrections, or even the full equations of motion in order to correctly reproduce the dynamics. A natural question is 
then what happens to the two metrics after non-linearity is reached: the evolution would then need to be calculated 
self-consistently in a non-linear theory for bigravity, which is beyond the scope of this paper. This problem is not new 
to Dark Energy models. There are cases such as growing neutrino cosmologies m in which the effect of non-linearities 
becomes important and needs to be taken into account also when dealing with the CMB predictions (55]. In that 
scenario, the fast growth actually leads to stable non-linear structures (which are a way to test the model rather than 
an argument to exclude it based on linear theory). 

Following the idea of Ref. (48], we then stop the evolution of perturbations at some cut off amplitude value h cut ss 1, 
when non-linearity is approximately reached. This prescription is applied to both h g and hf and has the effect to 
partially stop the ‘dragging’ of the second metric hf over the standard one h g . Such assumption is adopted here for 
simplicity, as for now our interest is to give a consistent estimate of how big is the impact of the growing mode on 
CMB spectra and tensor perturbations when non-linearity is reached. 
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In practice, since the growth of hf is very rapid in the early Universe, this is equivalent to fixing hf = h cut = 1 from 
the very beginning, with the consequence that the perturbations of the reference metric are not dynamical anymore. 
Nevertheless they still affect the tensor modes of the physical metric due to the coupling (321. The overall evolution 


as a function of the scale factor for two values of k and two different choices of the initial conditions is plotted in Fig(2] 
for both h g and hf , once the bound has been applied. Figure [2] also shows a similar behavior between a model in 
which hf starts saturated and one in which the cutoff value is reached during the evolution. It is also shown how, even 
though the cutoff is applied to both metrics, h g never reaches it during its evolution. Note that the full non-linear 
dynamics might produce other effects. For example, in the limit of scales smaller than the horizon, one finds the 
oscillating behavior in e-folding time h g ,hf ~ e lmN with eigenfrequencies: 


m = ±- 


Ur 


(46) 


This oscillating behavior is then always present in the sub-horizon solutions, overimposed to an amplitude modulation, 
as shown for the smallest scale in Fig. [2] (top right panel). In this case, setting hf to the non-linear cut-off value leads 
to a growing behavior plus a damping of the initial oscillations. On the other hand, the model with fine tuned initial 
conditions displays the oscillatory behavior expected from linear theory (461, and in this case the non-linear value 
hf ~ 1 is not reached, at least for Fourier modes corresponding to small scales. In this case the negative friction of 
the reference metric gets compensated by the positive friction from the physical metric. Nonetheless we expect our 
method to give a correct quantitative estimate of the observable effects of the growing mode on the CMB. 


IV. COSMIC MICROWAVE BACKGROUND ANISOTROPIES IN BIGRAVITY 


The Cosmic Microwave Background (CMB) was shown to be a powerful probe to test not only early time cosmology 
but also Dark Energy and Modified Gravity models |49| . In particular, in this paper we are interested in the effect that 
tensor perturbations in bigravity have on the CMB power spectra. At recombination, when photons are not anymore 
tightly coupled to baryons but decoupling has not occurred yet, electrons can be scattered simultaneously by photons 
coming from cold and hot spots. In presence of a quadrupole temperature anisotropy, the scattered photons will be 
linearly polarized and the CMB radiation will be characterized not only by its intensity, but also by its polarization. 
CMB polarization can be expressed in a tensor normal basis in Fourier space, in terms of E and B modes. While 
scalar perturbations can only produce an E mode (primordial) polarization pattern, tensor perturbations can feed 
both E and B primordial modes. Therefore any change in the evolution of tensor perturbations predicted in bigravity 
will affect the polarization spectra. 

At later times, polarization and temperature anisotropies are further modified during reionization. Reionization 
occurs at a much lower redshift, when the universe becomes partially ionized due to the formation of the first stars, 
allowing CMB photons to partially rescatter. The recombination and reionization eras correspond to peaks in the 
visibility function shown in figure [2] The visibility function g(t) = exp(— n)k (where k is the optical depth and k is 
its derivative with respect to conformal time r) gives the probability that a photon last scattered in the conformal 
time interval [r, r + dr]. Due to the importance of the coupling at relatively low redshift, the most important effects 
of tensor modes on the CMB are imprinted during the reionization epoch. 

In the following, we have only modified tensor perturbations, assuming that the contribution of scalar perturbation 
is small enough to be neglected, as scalar modes affect B mode polarization only indirectly, via lensing of E modes, 
at scales l > 150. Of course, if polarization is large enough, it might also feed back the scalar spectra. However, 
this seems a good enough first approximation to test the specific effect of the growing mode on the BB spectra. We 
implemented the tensor evolution equations in two publicly available Boltzmann codes, CAMB [50j and CLASS [5lj . 
and compared the results obtained in the various cases to verify their mutual consistency. 

As discussed in the previous section, our aim is to check the effect on the CMB spectra consistently, i.e. taking into 
account that, by definition, we cannot trust any result derived assuming linear perturbation theory when perturbations 
become non-linear. We then fix hf = h cut = 1 from the beginning and evolve only the ^-metric tensor h g , for which 
we will assume standard initial conditions: 

/ i, \ nT 

At = r T / S A s y — J , (47) 


with a fiducial tensor-to-scalar ratio tt/s = 0.05, a scalar amplitude A s = 2.21 x 10~ 9 and a tensor spectral index 
given by the self-consistency condition of single field sl ow-ro ll inflation ut = —(2 — tt/s !8 — n s)rT/ s/&, where the 
fiducial scalar spectral index is n s = 0.9645 (in section V A we explore the effects of changing the IC on the tensor 
modes.). For bigravity we choose the best fit model 0i = 0.48, 04 = 0.94 with £l c d m = 0.13, = 0.05, while for 
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FIG. 3: Effects of the growing mode in the /-tensors on the CMB. All the plots assume IBB with /3i = 0.48, /?4 = 0.94, 
fi m = 0.18 with flj, = 0.05, h = 0.67, Tt/s = 0.05 and the spectral index determined by inflationary self-consistency conditions. 
The tensor perturbations have been assumed to start at non-linear cutoff value hf = 1. 


the reference ACDM model we choose Planck 2015 TT, EE, TE+lowP marginalized values f52j, i.e. fl c dm = 0.26, 
fib = 0.05. In both bigravity and ACDM cases the fiducial optical depth is r = 0.079, corresponding to z re i 0 = 10. In 
order to test the IBB model, we compare the achieved spectra with up to date CMB observations, using Planck 2015 
data [53] for TT, TE and EE spectra, while for the BB spectrum we rely on WMAP [543 and BICEP1 [55] together 
with the joint BICEP2, Keck, Planck analysis [5fy. Figure [3] shows the tensor contribution to the CMB temperature, 
polarization and cross spectra for the fiducial bigravity model described above. Large angular scales are the most 
sensitive to the growing modes in the /-metric, because the coupling to the physical metric is only important at low 
redshifts, after recombination. On these scales the tensor perturbations give contributions to the power spectra that 
are too large to be compatible with CMB data. Even for the T and E polarization spectra, the tensor contribution 
in IBB bigravity overshoots the observed values by several orders of magnitude for £ < 60. Since both the scalar and 
tensor contributions to TT and EE spectra are positive definite, it is impossible that a reduced scalar contribution 
compensates the (large) tensor part in order to fit the data. The conclusion is that the cut-off on the growing mode 
is not enough to save the model. 
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FIG. 4: Non-solutions to the problem of growing modes in the reference metric: varying the IBB parameters across the 
degeneracy line (top left panel), the reionization redshift (top right) and the initial conditions for tensor perturbations in the 
physical metric (bottom panels). The remaining model details are the same as in Fig. [3] In particular, the reference metric 
perturbations are set initially to the non-linear cutoff value. All solid lines correspond to the standard values described in 
section m 


V. POSSIBLE SOLUTIONS TO THE GROWING MODE PROBLEM 

A first attempt to overcome the growing mode problem is to verify how much the effect on the spectra depends on 
the choice of the fiducial model. We investigated the effect of changing the bigravity parameters j3\ and choosing 
them close to the degeneracy curve (f24|) , the redshift of reionization z re io, the tensor to scalar ratio r^/s and the 
tensor spectral index ny characterizing the shape of the primordial power spectrum for tensor perturbations. These 
effects are shown in Fig. [4j showing that simple variations of these cosmological parameters do not offer a sufficient 
improvement. 

The observable impact on the CMB is produced in the reionization era because the coupling between the physical 
and reference metric’s tensor perturbations is only relevant at low redshift. Shifting 2 re i 0 in the range (5,15) has only 
a small effect on the BB spectrum. Further changes would spoil the predictions for EE and TE spectra and enter 
in contradiction with the Gunn-Peterson limit 152 ■ Higher reionization redshifts reduce the tension slightly, mainly 
because the perturbations in h g are smaller at earlier times and on smaller scales. Another attempt that does not 
work is to modify the initial conditions for perturbations of the physical metric. Varying the spectral index and the 
tensor to scalar ratio only has an impact on relatively high multipoles l > 30, on which the tensor perturbations are 
predominantly imprinted during recombination. For lower multipoles the evolution of h g is dominated by the coupling 
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to hf, which overshoots the initial conditions on large scales due to its large value. 

Concordance with observations can neither be achieved by varying the IBB para meters fa, fa. The fast growth of 
the hf tensor perturbations leads to growing h g tensors due to the coupling in (261 which is proportional to far fH . 2 . 
One might expect that a change in the betas could lead to higher values of r and, thus, a suppression of the coupling at 
early times (note that Ti 2 oc r 2 in RDE). Higher values of r are possible when lowering fa since ro = (1 — H tot o)//3i, 
where ro denotes the present value of r and one has r > rg during the entire evo lution. However, in order to fit 
observations we have to choose parameters being close to the degeneracy curve (241, i.e. fa/fa — const. As already 
mentioned, at early times the coupling term is then proportional to q g « 1/far (see Eq. 331. A smaller coupling 
would of course help in delaying or reducing the effect of the hf growing mode on h g . Since the evolution for large r 
is nearly independent of the coupling parameters (see Eq. 181, and since fa < 0.5, the ratio l/{far) cannot decrease 
much below the reference case. Figure [4] shows clearly that the effect of varying fas is negligible, at least when the 
coefficients remain along the degeneracy line. 

One loophole in the line of argument above is to leave the observational degeneracy curve ( |2~lj ) and add an explicit 
cosmological constant fa. This transforms the IBB model into a form of ACDM plus a small admixture of massive 
bigravity. The coupling coefficient q g does not explicitly depend on fa and reduces to: 


% = 


3 far 2 


fa + far 3 ' 


From the /^-Friedmann eq. (161 we find that r scales like fa 1 ^ 2 for small values of fa which leads to: 

9fa 


% - 


fa fa + 3fa3fa 


(48) 


(49) 


Thus, one is able to get an arbitrarily small coupling when choosing values of fa that are sufficiently smaller than fa. 
In this regime, the massless graviton dominates over the massive one and the cosmological evolution tends to that of 
ACDM. It is clear that such a model is not particularly interesting from a cosmological point of view because hardly 
distinguishable from the standard model and therefore here we will not investigate it further. 

In order to render the model viable without adding a cosmological constant, it is necessary to adopt a more radical 
solution. In the following subsections we explore how IBB can be reconciled with observations by fine-tuning the 
initial conditions in the reference metric perturbations, lowering the non-linear cut-off or modifying the theory. 


A. Changing the initial conditions for hf 


Our next attempt consists in checking whether fine tuning the initial conditions can compensate the effect of the 
growing mode on the CMB spectra, as illustrated in Fig. ([5]). We specify the initial conditions in terms of the growing 
solution in the radiation era 


^/(d) — ^/(in) (u/fli n ) ; hf — 3/ly(fl) . 


(50) 


found in section III B 
in 10 26 at Zi n = 10 iu 


We find that the initial conditions have to be fine tuned to zero to at least the level of one part 


in order to fit current limits on the BB spectrum. This choice of the initial epoch corresponds to 
an era before Big Bang Nucleosynthesis, which as already mentioned is a hard lower bound for the end of inflation. 
One can easily relate to earlier times in order to specify the IC at the reheating epoch, when inflation ends and tensor 
perturbations start growing. Table |T] extrapolates the result to the range of energies in which inflation might have 
ended. 


1 One can in general fi x hf and h/ independently for each wavenumber, but we restrict to the simpler choice ijoOli here. If more general 
IC are considered, a necessary condition for the growing mode to be sufficiently suppressed is that the time oerivative is small. This 
condition is sufficient as long as hf ( in ) is well below the cutoff value (see next subsection). 
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FIG. 5: BB (left panel) and TT (right panel) spectra using fine-tuned initial conditio ns a s described in section V A The initial 
amplitude and time derivative have been specified according to the growing solution (501 at a fiducial scale factor a ln = 10“ 10 . 
Bigravity BB spectra (red lines) only contain the primordial tensor contribution, while ACDM spectra (black solid lines) 
includes the contribution of both scalar and tensor perturbations. The evolution of tensor perturbations has been stopped 
whenever h g j = 1 is reached (cf. section HID I. 


_ BBN _ fiducial 1 GeV _ GUT 

Bounds on h f(in) < 10“ 19 < 10“ 25 < 10“ 31 < 10“ 82 

Scale factor a ~ 10“ 8 ~ 10“ 10 ~ 10“ 12 ~ 10“ 29 

Temperature T ~ O.lMeV ~ lOMeV lGeV ~ 10 16 Gev 


TABLE I: Upper bounds on hf extrapolated to different epochs using the growing mode, eq. (501. The limits shown are based 
on the results for BB spectra, which only depend on the tensor sector at low multipoles where the theory enters in tension with 
the data. Considering TT spectra would tighten these bounds, as can be inferred from the left panel of figure [5] 


As already noticed, the only way to generate naturally a very low level of tensor modes compatible with CMB 
without any cut-off is to assume that inflation ends at an energy scale not larger than 1 GeV. If however one assumes 
the non-linear cut-off, then the inflation energy scale bound can be relaxed. From Table |T] we see empirically that 
the initial condition for hf is related to the end of inflation energy scale T e (expressed in GeV) as hf( in \ sa 10“ 31 T“ 3 . 
Then one has 


h f (in) 



lGeV \ 

T e ) ’ 


(51) 


or T e k, 25 GeV, a more realistic scale range for low-energy inflation. 

A sufficiently small value of h'f might also be provided by a more exotic inflationary mechanism. During inflation 
no growing modes occur on the perturbations of the / metric. Some solutions, such as increasing the mass of the 
graviton at very early times, might naturally generate the low values needed to reconcile the model with observations 
(the problem of growing classical perturbations is common to ekpyrotic scenarios alternative to inflation [58j). Outside 
of these rather unconventional, although not impossible, cases, the conclusion we draw is that only very fine tuned 
initial conditions allow to reconcile bigravity with CMB observations. 

Even if inflation ends at a sufficiently low scale or a mechanism to suppress h'j^ exists, it has been argued by Cusin 


et al. |35J that non-linear corrections would spoil the small value of h'f. Although we will not investigate this issue 
further here, we note that the nature of the theory might protect the tensor modes against such terms. This is precisely 
what happens in the linear equation (261, in which the source term is highly suppressed. If a similar suppression 
occurs also on the non-linear source terms, the fine tuned initial conditions can render the model compatible with 
CMB observations (assuming that there are no additional complications in the scalar sector). 
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FIG. 6: BB (left panel) and TT (right panel) using a varying cutoff h cwt . The evolution of hf and h g is frozen when those 
reach, respectively, the cutoff value h cut . Bigravity BB spectra (red lines) only contain the primordial tensor contribution, 
while ACDM spectra (black solid lines) and observational data points contains the contribution of both scalar and tensor 
perturbation. 


B. Lowering the cut-off 

Another possible solution to mitigate the effect of IBB on CMB spectra is to assume that non linear effects begin 
to be not negligible before hf reaches unity; this can be thought as an effective way to treat the impact of non linear 
effects which, even if dominant when hf is above unity, can start to affect the evolution of the perturbations even for 
lower values. Therefore, in Figure [b] we show the behavior of TT and BB spectra for different values of h cut at which 
we freeze the evolution of the metric perturbations. One would need to suppress the cutoff scale by at least three 
orders of magnitude to reconcile theoretical BB spectra with currently available limits and possibly an even lower 
value to make the TT spectrum acceptable for current data, once the scalar contribution is taken into account. 

It seems contrived that non-linear effects might play a role at such small values of the cutoff. However, theories 
of massive gravity are known for having strong non-linear effects in certain limits, such as the Vainshtein mechanism 
150] , As long as numerical results for the non-linear evolution of tensor perturbations are not available, we must 
contemplate the possibility that the cut-off could be lower than unity and even significantly lower. 


C. Modifying the theory 

The growing modes in the reference metric could possibly be reconciled with CMB data by a suitable modification 
of the theory. Here we explore a phenomenological modification in which a redshift dependence of the parameters is 
assumed; this kind of behavior might be achieved in generalized massive gravity models ED, where a time dependence 
of the mass parameters is introduced without the addition of any new dynamical degree of freedom. Our modification 
consists on setting /3\ = P 4 = 0 until a certain switch redshift z s is reached. Then bigravity becomes active and the 
evolution described in Section |III| is switched on. The considered background evolution instead is the one produced 
by bygravity at all redshift, as IBB well approximates the standard background, which should take place at z > z s , 
at early times. 

The results in Figure [ 7 ] show how switching on bigravity at approximately the redshift of matter-radiation equality 
(z s ~ 10 3 ) can produce an acceptable BB spectra when comparing with current data. In the TT case instead, the 
contribution of scalar modes to the spectrum can possibly lead to the necessity of an even lower value of z s . 

There are additional ways in which the theory might be modified while retaining the original field content of 
bigravity. In the following we will describe these possibilities, although addressing them in detail will be left for future 
work. A possible modification is to allow for branches different than IBB. So far, only branches in which r evolves from 
r = 0 or r —> 00 were considered. In [2p. all remaining cases due to a non-viable behavior were excluded. However, 
some of those conditions based on expectations of a standard cosmological evolution, like an expansion at all times 
and the existence of a matter/radiation dominated era in the asymptotic past. It might be interesting to study these 
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FIG. 7: BB (left panel) and TT (right panel) using a varying initial redshift z a . If 2 > z s hf follows the same standard equation 
as h g , while for z < z s the evolution described in Section III is switched on. Bigravity BB spectra (red lines) only contain the 
primordial tensor contribution, while ACDM spectra (black solid lines) and observational data points contains the contribution 
of both scalar and tensor perturbation. 


disregarded branches. Additionally, both metrics are usually assumed to be FLRW at background level, even though 
this is assumed only for simplicity. In Ref. |62] the authors considered more general types of metrics which might 
lead to interesting evolutions at background and linear level and could also have an impact on the tensor evolution. 

An additional possibility is to modify the coupling to matter. Even though both metrics are a priori equally footed, 
we let only one of them couple to matter while the remaining metric stays unobservable. An additional coupling of 
the reference metric to matter would influence the tensor perturbation and might be able to tame the fast growth 
(see e.g. Refs. [63] !63 for further discussion on bi-metric couplings) . An additional coupling of the same matter 
Lagrangian to / M „ is not possible as it will generally reintroduce the BD ghost (65fl67| . Even though this will not 
happen if a different matter Lagrangian (an unobservable dark sector) is coupled to / M „, we will usually meet a new 
fine tuning [68] that would make the theory less appealing. One way out would be a coupling through a new composite 
metric that is constructed such that it avoids the BD ghost [65] ■ This choice would lead to viable, self-accelerating 
backgrounds [69] but still does not yield a realistic cosmological evolution at the linear level Emm. 


VI. CONCLUSIONS 

We have analyzed the behavior of tensor perturbations in the infinite-branch bigravity (IBB) model and the signa¬ 
tures they produce on the CMB. In this model the reference metric contracts at early times, causing tensor pertur¬ 
bations in this metric to grow rapidly. These modes have ample time to grow since the end of inflation. However, 
a consistent analysis of linear perturbations and CMB spectra cannot include the regime in which perturbations be¬ 
come non-linear. We then assume that this growth stops when perturbations become non-linear, with an amplitude 
saturated at a value of order unity. The coupling between the two metrics produces in turn a growth of the tensor 
perturbations in g. If the coupling is weak enough at early times, the growing mode will in principle propagate to the 
physical metric only after some time. Our first objective has been then to check whether this effect is late enough 
to keep the spectra of tensor perturbations compatible with present CMB data. Our conclusion is that, even when 
perturbations remain below order unity, they are still large enough to have a twofold impact on CMB spectra: first, 
the tensor modes provide a large contribution to the TT, EE and TE spectra which is orders of magnitude larger 
than the scalar contribution; second, the tensor modes induce a strong B-mode polarization on the CMB. Both effects 
dominate on the largest angular scales and are incompatible with observations from CMB experiments. 

Varying the IBB parameter (/3 1 , with (3^ being derived from it via Eq. ©) or other cosmological parameters offers 
little help in reconciling IBB bigravity theory with observations. We further explore five scenarios in which the theory 
might be rendered viable: 

1. Lowering the energy scale of inflation. If the energy scale of inflation is very small, around 1 GeV, the tensor 
modes are naturally suppressed and the growing mode has less time to grow until recombination or reionization: 
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the combined effect makes the IBB model acceptable without any change. If one invokes the freezing of the 
growing mode when it reaches non-linearity, then the inflationary energy scale can be increased up to 25 GeV 
roughly. Such a scale is much lower than the one predicted in simple slow roll inflation but could in principle 
be achieved in alternative scenarios [551 . 

2. Fine-tuning the initial conditions. If h'j is very small at early times, the perturbations will not have reached 
the non-linear value at late times, when the coupling to g becomes important. This requires a fabulous degree 
of fine-tuning, of one part in 10 26 at z = 10 10 (10 28 when TT modes are considered). For this solution to work 
beyond the linear approximation, it is necessary that non-linear sources in the equation for hf are suppressed 
at early times. Even if stable agains non-linear corrections, fine tuning the initial conditions for hf seems a 
highly ad-hoc requirement for the theory in the absence of a mechanism, perhaps generalizing inflation, able 
to naturally produce such small initial conditions. As mentioned in the previous point, this fine-tuning occurs 
naturally only in the case of very low-energy inflation. 

3. Lowering the non-linear cutoff. IBB becomes safe if hf stabilizes at a value smaller than unity due to non-linear 
effects. This requires the non-linear effects to act at most when hf < 10 -3 (and < 10 -4 when TT modes are 
considered). 

4. Adding a cosmological constant. In that case the theory still describes a massive graviton, although it would 
not be responsible for the acceleration of the universe. From a purely cosmological point of view, the model will 
be very similar to ACDM. 

5. Modifying the theory. One possibility is to allow for a time dependence to the theory parameters /3i, /?4 in the 
tensor perturbation equations. This phenomenological parametrization of modified IBB allows to satisfy CMB 
data if the parameters become non-zero only after z = 1000. 

While our assumption stops the growth of perturbations when they become non-linear, a full analysis would require 
an understanding of the actual non-linear behaviour of perturbations, which could be used to exclude or validate such 
scenarios in a fully consistent way. 

Finding a modification of bigravity that overcomes the difficulties of the growing modes could be possible within the 
framework of generalized massive gravity m ■ In this class of theories, the interaction terms are given an additional 
dependence in the Stiickelberg fields, which allows the couplings to vary over time without introducing additional 
degrees of freedom. Another possibility in this direction is allowing a composite coupling that involves both tensor, 
possibly on an equal footing. Theories with additional degrees of freedom, such as scalar-bitensor or multigravity, 
might as well prove useful to solve the problem of growing modes. More exotic modifications of the theory remain to 
be explored. 

Constructing a viable theory of massive gravity has proven to be a challenge. Only after eight decades could the 
linear Fierz-Pauli theory be generalized to a ghost-free non-linear completion, albeit one that forbids any interesting 
cosmological solution. This difficulty could be overcome by giving a kinetic term to the reference metric, allowing the 
existence of accelerating cosmologies at the price of two additional degrees of freedom, corresponding to a massless 
tensor. Of all the five-parameter set of bigravity theories, only the two-parameter IBB family is able to accelerate the 
universe with neither a cosmological constant nor scalar instabilities. Yet, such a theory is affected by growing modes 
that generically spoil the predictions of the cosmic microwave background. The results presented here represent a 
setback for the simple and appealing self-accelerating bigravity paradigm, a paradigm that, unless saved by non-linear 
effects or a tiny amplitude of the initial conditions, will have to be abandoned. 
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